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Abstract 

The addition of minute amounts of chemically inert polyacrylamide polymer to liquids results in 
large instabilities under steady electro-osmotic pumping through 2:1 constrictions, demonstrating 
that laminar flow conditions can be broken in electro-osmotic flow of viscoelastic material. By 
excluding shear and imposing symmetry we create a platform where only elongational viscoelastic 
instabilities, and diffusion, affect mixing. In contrast to earlier studies with significant shear that 
found up to orders of magnitude increase in mixing we find that inclusion of polymers excites 
large viscoelastic instabilities yet mixing is reduced relative to polymer-free liquids. The absolute 
decrease in mixing we find is consistent with the understanding that adding polymer increases 
viscosity while viscoelastic flows progress towards elastic turbulence, a type of mild (Batchelor) 
turbulence, and indicates that electro-osmotic pumped devices are an ideal platform for studying 
viscoelastic instabilities without supplementary factors. 
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I. INTRODUCTION 



This paper explores how mixing is affected by adding high molecular weight polymers 
to electro-osmotically pumped microffows. Adding low concentrations of polymer to fluids 
moderately increases viscosity and dramatically increases elasticity, allowing microffows to 
undergo viscoelastic instability. Based on the relatively small increase in viscosity, the 
striking qualitative change in stability, and early experiments that demonstrate up to orders 
of magnitude improvement in mixing there is an expectation of dramatic increase in mixing 
by incorporation of polymers. By excluding shear, present in all mechanically pumped 
cavity flows, we find that contrary to such expectations absolute mixing is reduced by 
adding polymer to induce viscoelastic instability. Our findings demonstrate that viscoelastic 
instabilities do not, in themselves, promote mixing and suggest that electro-osmotic flows 
are an ideal platform for isolating and studying behavior of viscoelastic instabilities and 
elastic turbulence. 

Viscoelastic flows are of practical interest for lab on a chip applications as many biological 
materials of interest for biomedical testing, such as saliva [I] and blood [2], are viscoelastic 
in nature as are many other suspensions [3] . Several devices sort particles using laminar flow 
interactions suggesting that viscoelastic induced instabilities can have potential unde- 
sirable effects as more concentrated and complex fluids are investigated; it has already been 
noted that viscoelasticity in saliva testing negatively effects performance |T]. In addition to 
the detrimental effects of instabilities there is suggestion that enhanced mixing may occur by 
intentionally spiking samples with polymers [5], thereby causing instability and subsequently 
increased mixing. As chemical reaction rate and completion depends on good mixing and 
as many clinical, environmental, pharmaceutical, and other diagnostic applications require 
rapid and complete reactions, in particular in the presence of small concentrations of sample, 
efficient mixing is an important aspect to device performance. 

In microfluidic devices the scales involved (typically 0(100) /xm) are such that diffusive 
processes are slow, yet inertial fluid instabilities are suppressed due to the small scales leading 
to small Reynolds numbers. This intermediate "slow scale" regime makes mixing a key issue 
in microfluidics [6\ and various methods have been investigated in order to enhance mixing, 
including using bubbles [7], electrohydrodynamic instability [8J, flow focusing and bas- 
relief patterning JTU]. These approaches generally rely on creating inhomogeneities to the 
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working fluid or the microchannels. Viscoelastic fluids can exhibit flow instabilities [TT] 
when polymer coils stretch and bend along stream lines which creates hoop-stress, allowing 
polymers to cross steam lines and thereby disrupting laminar flow [T2] . 

These viscoelastic instabilities progress to elastic turbulence [13], a recently established 
turbulence that arises from elastic effects that can be excited at arbitrarily low Reynolds 
numbers. Use of viscoelastic fluids to promote mixing has been demonstrated in both mi- 
croscaled Dean flows [14] , where chaotic flow was established and where earlier study of a 
macroscaled version of the serpentine flow geometry demonstrated true elastic turbulence 
[5]. The establishment of viscoelastic turbulence JT3] and efficient mixing [5] along with 
favorable scaling properties of viscoelastic instabilities, which become easier to excite with 
reducing dimensions, suggests that the use of viscoelastic fluids may be a promising route to 
promote mixing in microfluidic devices. Furthermore, it is well know that solution viscosity 
increases with polymer concentration and at the overlap concentration, where neighboring 
coils impinge on each other, viscosity is only double that of the polymer-free solution [13] . I n 
contrast to this relatively mild quantitative change in viscosity (and thus decrease in diffu- 
sion) elasticity effects have dramatic qualitative impact with polymer concentrations as low 
as ~ 10 ppm leading to both large drag reduction in high Reynolds number hydrodynamic 
turbulent flows [15] and sizable instability and elastic turbulence in low Reynolds number 
flows [5]. 

Several studies (see, for example, [HI [TBI HZ]) have found that viscoelastic instabilities 
can be excited in microdevices, however additional deformations of the flows - most notably 
shear - exists in prior work and highly viscous solvents are used, which is atypical for many 
lab on a chip applications. Here we use electro-osmotic flow (EOF) in order to exclude 
shear and solely induce elongational viscoelastic instabilities by driving low viscosity fluids 
through sudden 2:1 constrictions in linear microchannels; we compare mixing in unsteady 
viscoelastic flows relative to mixing in polymer-free solutions in order to measure mixing 
effects. 

Mechanically driven cavity flows always have a significant shear component to the flow 
due to no-slip conditions. In contrast EOF has effective slip boundary conditions and shear 
is confined to the Debye layer close to the channel boundary [I]; flows are shear-free external 
to the Debye layer as the electric fields that the flow follows are purely elongational [T5] . 
Further, for large polymers, entropic repulsion from the walls reduces interactions of the 
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FIG. 1: (a) Photograph of the microfluidic device, and (b) detail of the micro channels with mixing 
units. We electro-osmotically drive fluids from the top and bottom wells to the left side well. Fluid 
in the top well contains fluorescent dye, and mixing is optically observed by exciting the dye 
with an external 532 nm laser and collecting video data. Each 100:200 /urn constriction/expansion 
unit creates extensional flow which induces instabilities in viscoelastic fluids with high enough 
Deborah number. Due to the slip conditions of electro-osmotic flow and symmetry non- viscoelastic 
instability mixing mechanisms are minimized. We find flow rates are not affected by the presence 
of polymer in the dilute regime and drive all flows with a voltage of 200 V, resulting in equal flow 
rates and comparable flows. 

Debye shear regime with polymers |IH]. Our pumping scheme therefore excludes the shear 
component that is ubiquitous to mechanically pumped flows. In addition to deformations 
caused by shearing there are possible Lagrangian chaos effects [20], where flow fields result 
in chaotic trajectories which can lead to rapid mixing. We stack identical constriction units 
equally spaced down a microchannel, imposing symmetry as to reduce possible Lagrangian 
chaos effects [2"T] . 

II. EXPERIMENTAL SETUP 

See Ref. [22] were we demonstrate viscoelastic instabilitites under steady electro-osmotic 
pumping and discuss our device in more detail. Briefly, the experimental details are as 
follows: a microfluidic toolkit (/iTK) from MicraLyne Inc. [B] is used to impose EOF flow 
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in glass microchannels etched 20 fim deep. TAMRA dye, a fast diffusing dye with D m 270 
/xm 2 /s [21], fluorescently labels one, otherwise identical, stream of fluid allowing visualization 
of flows. A dyed and an un-dyed fluid stream are brought together in a cross-slot junction 
and flow downstream through 200:100 /im constrictions (see Fig. [I] for images of the device 
and microchannels). Fluid flow details are captured using a a Casio Exilim EX-F1 camera in 
video mode (30 fps), by attaching the camera to the /iTK via a custom machined attachment 
and a Scoptronix Maxview Plus optical adaptor (camera-microscope attachment). 

Our camera is not scientific grade and comparison of CCD intensity with true (incident) 
laser power reveals nonlinear characteristics; cubing the CCD values results in a linear 
response profile allowing the optical intensity as measured to be scaled and related to the 
dye concentration. From video data space-time (xt) diagrams are created by taking a vertical 
cut from each video frame across the channel, horizontally centered in the constriction or 
expansion. Due to chemical etching side walls exhibit curvature, varying the channel depth 
near the edges, and we remove 20 /zm nearest the channel edges before analysis. Runs are 
nominally 100 s in duration, but we exclude initial transient flow to prevent start up effects 
from influencing results. 

Polymer solutions were made by adding high molecular weight (18 x 10 6 Da) polyacry- 
lamide from Polysciences Inc. to a 20:80 vol.% methanol:water mixture; methanol is used 
to prevent aging during storage [25]. The overlap concentration c* is ~ 300 ppm, and we 
work in the dilute regime (c < c*). 

Flow velocities in EOF are given by v — HeofE, where /ieof is the electro-osmotic 
mobility and E the applied electric field. We find that below the overlap concentration 
the electro-osmotic mobility is not significantly modified from the pure solvent case, with 
/xeof ~ (5.6 ± 0.5) x 10 -4 cm 2 /Vs here, a result that is predicted by the entropic wall 
exclusion effect and studied in Ref. [19] . Our flow velocities are on the order of mm/s, and 
Re << 1, as is typical of standard EOF microflows. 

III. RESULTS AND DISCUSSION 

Addition of polymer allows elastic effects to induce instabilities for sufficiently high Deb- 
orah number De above the instability threshold (> 1/2 for elongational flows [26]), the 
dimensionless number that characterizes viscoleastic flows. The Deborah number is given 
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by 

\V XzimmP-EOFE . . 

De=- « . (1) 

-^Char -^Char 

We estimate the relaxation time A of our liquid using the Zimm estimate Xzimm ~ 0.03 s 
and Lchar ~ 12 /mi for our channels; see Ref. |22j. For our setup and polymer De > 1/2 for 
driving voltages > 20 V, and here we operate well above the threshold at 200 V to ensure 
highly unstable flows. As we find //eof is not modified by the presence of polymer in the 
dilute concentration regime all flows are comparable. 

The strategy to promote instabilities discussed here is simple. The addition of polymer 
is facilitated if pre-solvated polymer mixtures are added to working fluids. The use of high 
viscosity solvents or concentrated polymer is not required: low viscosity, lightly polymer 
doped liquids undergo instability under typical electro-osmotic flow conditions. 

The mixing efficiency is found by measuring the homogeneity of a fluid domain [27], and 
here is calculated using the first moment mixing index 

Ml = (|C ' 5|) (2) 
c 

where c is the dye concentration and c the average concentration; for unmixed streams 
Ml = 1 and mixing is perfect for Ml = (homogenous). 

In "normal" microfluidics smooth laminar flows are observed, as seen in Fig. [2] for our 
polymer-free solutions. As the fluid flows downstream diffusion across the channel will 
eventually result in homogenously mixed fluid (not show). We observe that mixing values for 
polymer-free flows scale approximately exponentially as a function of distance downstream; 
while perfect exponential scaling is not expected calculating Ml using the theoretical erf- 
profile [28] for laminar mixing in uniform channels results in an approximate exponential 
drop off [45j. 

Addition of polymer results in large instabilities, as can be seen in Fig. [3] and Fig. |4} 
The fluctuations move material transverse to the flow direction, spanning almost the entire 
channel. There is no dominant frequency in the flow, indicating chaotic flow [Ti] . and the 
overall pattern displayed are striations with "comb-like" teeth randomly spaced. Increasing 
the concentration from 32 to 64 ppm results in larger fluctuations and tooth spacing. As 
the solutions move downstream through additional elongation regions additional striations 
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FIG. 2: Laminar flow in the first constriction. With no polymer added creeping laminar flow with 
steady diffusion is observed (a). A space time diagram (a) is constructed by taking a cross channel 
intensity profile across the middle of the constriction and is 60 p by 75 s. Note that 20 fim is 
trimmed from the edges (where wet etch rounding may lead to aberrations), and that transitory 
start-up flow is removed, from all xt-diagrams here. The cross channel intensity profile 0.6 s from 
the start of the xt diagram is shown in (b). An erf- like profile is observed in (b), as expected in 
coflowing two-stream mixing. 

are added, as seen in Fig. [5j however the general pattern persists with diffusion appearing 
to be the dominant mixing mechanism. 

Absolute mixing is reduced by reduced diffusion, and increased by instability induced 
deformation of the liquid. There is therefore a tradeoff involved in adding polymer to the 
solution which will allow viscoelastic instabilities yet reduces diffusion. By measuring mixing 
(Ml) in the second expansion we find that mixing is reduced as concentration is increased (see 
Fig. [6]), demonstrating that the effect of instabilities in promoting mixing is not dominant 
for mixing our dye. As our dye is a small molecule (~ 1.6 nm hydrodynamic diameter [24J) 
our mixing experiments is a stringent - but realistic - test. Further, as we trim the edges 
of the xt-diagrams to reduce image processing artifacts we possibly artificially increase the 
first moment mixing index Ml calculated for unstable flows, as material with high/low dye 
concentration can be injected from these regions, increasing the variance seen. 

Note that in contrast to high Reynolds number turbulent flows an enhanced diffusion 
rate, indicative of broad-band excitation of instability scales, is not observed |29j supporting 
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FIG. 3: Photomicrograph of unstable flow in the second expansion for a concentration of 64 ppm 
high molecular weight (18 x 10 6 Da) polyacrylamide polymer added. It is apparent that large 
instabilities are excited that display a dominant scale set by the microchannel. Dashed outlines 
indicate field of view and microchannel outline; the 160 (im long solid line indicates where the 
space-time diagram slice was taken for the expansion regions. The instabilities in the constriction 
regions can also be seen here at the sides of the image; regions are centered for data collecting. 

the limited bandwidth visual evidence (see Fig. [3]) that large scale instabilities dominate the 
instability spectrum. The lack of enhanced diffusion confirms that small scale instabilities 
below the camera resolution are not excited to a significant extent. 

Finding the first moment mixing index, Ml, for concentrations of ppm, 32 ppm, and 64 
ppm as the fluid moves downstream reveals that absolute mixing is reduced in viscoelastic 
solutions relative to the polymer-free sample; Fig. [7j We also observe an increase in the rate 
of change in mixing downstream for viscoelastic solutions. By looking at the flow at different 
constriction units (N=2, 4, 8) the mixing length is calculated and we find a mild decrease 
with raising polymer concentration; the mixing length reduces from 4.5 for ppm to 4.2 
for 32 ppm and finally 3.2 for 64 ppm. For diffusion from an ideal sharp 2D interface [28] 
changing viscosity simply shifts calculated Ml by a constant for our parameters, indicating 
that the change in mixing length can be attributed to the action of instabilities. 

For an uniform microchannel we calculate Ml values roughly twice as large as we 
observe for polymer-free solutions in our channels corrugated with constriction/expansions. 
This enhanced mixing has been previously observed for electro-osmotic flows through con- 
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FIG. 4: Space-time diagrams of unstable flow in the second constriction. Adding 32 ppm (a) and 
64 ppm (b) polymer results in instabilities as fluid flows through the constriction/expansion units. 
Visual comparison indicates that instabilities more rapidly fluctuate across the channel for lower 
concentration, and that instabilities move more material across the channel with increased polymer 
concentration. The diagrams are 60 /im by 80 s, and have been scaled here to increase the contrast. 

strictions/expansion units where » 2x enhancement was found [30J. This enhancement is 
an open problem currently being investigated [31], but is thought to be due to the modula- 
tion of width enhancing diffusion [3T| due to convective cross flow movement (e.g. periodic 
movement perpendicular to the main flow direction) [46J. 

We note that Lam et al. Ref. [32J create instabilities by driving viscoelastic fluids 
through abrupt 8:1 constrictions using a syringe pump, and therefore have a mechanically 
driven device similar to our device. They observe enhancement of mixing in their liquids 
over deionized water at equal flow rates, demonstrating that absolute mixing increase can be 
obtained over low viscosity Newtonian liquids. However, they use dis-similiar fluids of differ- 
ent viscosity, shear-thinning, and elastic properties making interpreting their results difficult 
as flow focusing, shear-thinning [33], and viscous thread folding [31] or other coflowing ef- 
fects are physical mechanisms that all may play important roles in addition to viscoelastic 
instabilities in promoting the mixing they observe. 

We have shown that viscoelastic instabilities alone will not lead to rapid mixing, in 
contrast to earlier studies [5j [HI HZ] that found up to orders of magnitude improvement 
in mixing for mechanically driven channel flow. However, significant shear stretching is 
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FIG. 5: Mixing evolution. From the top to the bottom 64 ppm sample passes from the 2 nd (a) 
to the 4 th (b) to the 8 th (c) constriction. The effects of diffusion are apparent, with a reduction of 
contrast between the 2 nd (a) and the 4 th constriction (b), and finally to a visually homogeneous 
mixture by the 8*^ constriction (c). In addition it can be seen that strong cross channel fluctuations 
exist, but the pattern is highly similar and a comb- like pattern persists between the 2 nd (a) and the 
4 th (b) constriction, with addition of "teeth" but without significant pattern changing deformation. 
This persistence is characteristic of chaotic flows undergoing symmetric oscillatory forcing. The 
diagrams are 60 /im by 60 s and have not been scaled, allowing comparison. 

present in most prior work due to no-slip conditions. In a study of electrokinetic driven 2D 
sheet flow [35] , which also excludes no-slip induced shear, viscoelasticity was found to mildly 
suppress mixing for flows below the instability threshold, a finding that is consistent with 
other studies that find either mild increase |36j or decrease [37] in this regime where polymers 
are only weakly distorted. Additionally, in Ref. [35J they show that above the instability 
threshold stretching fields are reduced in viscoelastic fluids. By comparing a Newtonian and 
a highly shear-thinning/viscoelastic fluid they find enhanced stretching and mixing for the 
shear-thinning/viscoelastic fluid, indicating that mixing trends are correlated to stretching 
field magnitude. Our finding of reduced mixing, combined with the reduced stretching 
fields observed in Ref. [35], shows the correlation between stretching fields and mixing 
holds for viscoelastic fluids above the instability threshold. Our results which demonstrate 
reduced mixing in unstable viscoelastic flows are consistent with the understanding that 
viscoelastic fluids driven well above the viscoelastic instability threshold progress towards 
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FIG. 6: Mixing as a function of concentration in the 2 nd expansion unit, as measured by the first 
moment mixing index Ml (Ml — > for perfect mixing and Ml —> 1 for no mixing). Using space- 
time diagrams mixing is measured for c = 16, 32, 64, 128, and 256 ppm (the overlap concentration 
is c* ~ 300); error in Ml is estimated as < 5%. The mixing index slowly raises with concentration 
in the dilute regime; this indicates that while instabilities are promoted by the presence of polymers 
absolute mixing is reduced (diffusion is hindered by the presence of polymers). The lack of enhanced 
diffusion demonstrates that small scale instabilities are not excited to any significant extent. 
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FIG. 7: Mixing as a function of constriction unit, as measured by the first moment mixing index 
Ml. Using space-time diagrams mixing was measured downstream of the initial junction at N=2, 4, 
8. The mixing index scales exponentially, and the mixing lengths mildly reduce with concentration 
(4.5 for ppm, 4.2 for 32 ppm, and 3.2 for 64 ppm). The vertical shift for polymer doped fluids 
relative to undoped fluid is due to reduced diffusion limiting the mixing, as increased polymer 
concentration reduces the diffusion constant by increasing viscosity. 

elastic turbulence which has been shown to be a realization of Batchelor turbulence [38], [39] . 

In Batchelor turbulence it has been established that the velocity spectrum decays rapidly 
(the wavenumber k scales as k~ x , where x has experimentally been found to be in the 



0.8 
Jl 0.6 




0.2- 



o.o4- 
o 



43 0.4- 



* 

o 



0.6 



11 



range 3.3 to 3.6 |39j) which leads to weak coupling between scales, in sharp contrast to high 
Reynolds number turbulence, and the resulting flows are random in time yet spatially smooth 
and strongly correlated on the integral scale [3S]. This results in instabilities dominated by 
large scale fluctuations, which are set by the system size [39], and therefore there will remain 
a mismatch between the smallest scale one wants homogeneity on and the instability which is 
intended to drive mixing. Note that the smooth flow characteristics of Batchelor turbulence 
are consistent with our observations (see, for example, Fig. [3]) and, without some other 
mechanism (such as shear) to further modify flow, are at odds with enhanced mixing. 

To create efficient mixing a bridge is needed between the relatively large scale set by 
viscoelastic instabilities and the smallest scales over which homogeneity is required; diffusion, 
Lagrangian chaos, or deformations such as shear can bridge these scales. In the "knead and 
fold" metaphor of mixing [40J viscoelastic instabilities can provide the fold operation that 
creates large deformations of material to create striations of dis-similiar material, however 
additional deformations are needed to thin striations allowing diffusion to finally homogenize 
the fluids. 



IV. CONCLUSION 



In summary, the addition of high molecular weight polyacrylamide to liquids leads to 
striking viscoelastic instabilities in extensional electro-osmotic flow thereby upsetting the 
laminar flow typical of micro channels. The breaking of laminar flow for dilute low viscosity 
viscoelastic fluids has implications for device design as laminar flow is often assumed in 
microflows; as lab on a chip applications move towards more complex fluids this assumption 
must be checked. Our results show that for flows of lightly polymer-doped low viscosity 
solutions through unoptimized geometries mixing is reduced relative to polymer-free solu- 
tion, despite the excitation of dramatic large scale viscoelastic instabilities. This reduction 
is in contrast to earlier studies that found up to orders of magnitude increase in mixing, 
demonstrating that the role of shear or other deformations of the flow is crucial in promoting 
mixing. Electro-osmotic flows can be used to isolate viscoelastic instabilities, and combined 
pressure/EOF driven flows [JT] allow controlled application of shear enabling detailed in- 
vestigation of flow evolution and mixing. In addition to fundamental studies possibility for 
application exists for the observed highly unstable, yet diffusion limited, flows in applica- 
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tions which depend on large interfacial area; for example, in membraneless fuel cells jl2] or 
for use in combination with the liquid junction potential |43j. 
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The linear norm of residuals is < 0.01 over the range of mixing units tested, indicating expo- 
nential scaling is a close approximation of mixing behavior here. 

Naively considering the sheaves of flow lines suggests that the total distance an average particle 
travels as it moves downstream in a corrugated channel will be larger than the minimum 
distance (e.g. channel length); fluid off center will effectively "see" a longer channel and thus 
experience more diffusion, and one should expect enhanced overall diffusion in a given length 
relative to the centerline flow. 
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